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I review some open problems on the ever-growing field of non-equilibrium phase transitions, pay- 
ing special attention to the formulation of such problems in terms of Langevin equations or, equiv- 
alently, field-theoretical descriptions, and their solution using renormalization group techniques. 
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I. INTRODUCTION 

The introduction of the Renormalization Group (RG) 
ideas and their application to the study of critical phe- 
nomena constitutes one of the milestones in the spectac- 
ular development of the Statistical Physics during the 
last quarter of the twentieth century. The RG proved 
to be not only a powerful analytical tool to deal with 
complex problems {i.e. , problems with many different 
relevant scales) , but also a conceptually beautiful and ele- 
gant theory, with a huge range of applicability. While RG 
ideas applied in (discrete) real space have helped to shed 
light on many problems |l| , it has been in the framework 
of (continuous) field theoretical descriptions where, com- 
bined with perturbative methods, the RG has reached its 
most successful expressions. In particular, given a con- 
tinuous field theoretical representation of a given statis- 
tical system at thcrmodynamical equilibrium, the iden- 
tification of its critical points with fixed points of a con- 
veniently defined RG transformation permits to obtain 
(perturbative) analytical expressions for the correspond- 
ing critical exponents. At the same time, by allowing 
to distinguish relevant from irrelevant ingredients in a 
rather systematic way, the combined use of field theories 
and the RG, has permitted to obtain elegant classifica- 
tions of equilibrium critical phenomena and put under 
firm basis the concept of universality [|[^. The repre- 
sentation of the Ising model universality class by the Xcj)^ 
theory and its subsequent renormalization constitutes a 
paradigmatical instance 

Given this rather satisfactory scenario, theoreticians 
started wondering whether also critical phenomena oc- 
curring in systems away from equilibrium could be at- 
tacked using similar tools. In the lack of a well es- 
tablished theory for non-equilibrium phenomena, it is 
not straightforward to extend the equilibrium formalism 
to deal with non-equilibrium problems (for example, in 
these cases there is no partition function to be renormal- 
ized) . The best way to do so turned out to be expressing 
such problems in terms of Langevin equations describing 
the underlying dynamics at a (continuous) coarse grained 
level. This procedure is valid not only to study general 
non-equilibrium processes but also relaxation to equilib- 
rium states (the model A and B for the non-conserved 
and conserved relaxation dynamics of the Ising class are 
prototypical examples). In some cases Langevin equa- 
tions representing given microscopic processes can be 
derived analytically using different techniques (among 



them: Fock space formalism combined with coherent 
state representations jsj, and Poisson transformations 
[Q), while in many other cases they are just postulated 
from a phenomenological ground, by respecting what are 
considered a priori to be main symmetries, conservation 
laws, and other relevant dynamical constraints. 

Experience teaches us that the richness and variety of 
phase transitions occurring away from equilibrium is by 
far much larger than that of equilibrium, and that in 
many cases it is very difficult to decide a priori what 
are the essential ingredients to be put into a sound 
Langevin description. Therefore, developing rigorous and 
systematic techniques envisaged to derive coarse grained 
Langevin equations from microscopic models is a high 
priority task within this context. On the lack of such 
general approaches one has to rely many times on phe- 
nomenological approaches. 

Before proceeding, let us remark that any arbitrary 
Langevin equation can be written as an equivalent 
Fokker-Planck equation [^IJ], and its solution expressed 
in terms of a generating functional (or equivalently and 
effective action) obtained as a path integral represen- 
tation of the stochastic Langevin process. Therefore, 
in what follows "Langevin equations" , "Fokker-Planck 
equations" , "generating functionals" , or "field theoretical 
actions" (Liouville operator) will be used interchangeably 
(see [1101,1). 

In the forthcoming sections we report on a variety of 
interesting non-equilibrium systems, and present a list of 
open problems within this field. 



II. THE DRIVEN LATTICE GAS (DLG) 

The DLG is a variation of the kinetic Ising model with 
conserved dynamics, in which transitions in the direc- 
tion (against the direction) of an externally applied field, 
E, are favored (unfavored) [^|j9|Jl0[|, while transitions per- 
pendicular to the field are unaffected by it. The external 
field induces two main non-equilibrium effects: (i) the 
presence of a net current of particles along its direction, 
and (ii) strong anisotropy. At high temperatures, the sys- 
tem is in a disordered phase while lowering the temper- 
ature there is (for half-filled lattices) a continuous tran- 
sition into an ordered phase with high and low density 
aligned- with-the-field stripes. Elucidating the DLG crit- 
ical properties is an important issue in the way to ratio- 
nalize the behavior of non-equilibrium systems. The fol- 
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lowing Langevin equation was proposed some years back 
based on phenomenological arguments Illl.nGl : 



+r||V|0- aV 



A. 

6 



(1) 



where (j) is the coarse grained field, 77 is a conserved Gaus- 
sian noise, and where the cubic term (a dangerously irrel- 
evant variable is kept in order to ensure stability 0. 
The fact that some of the predictions derived analyti- 
cally from this equation (for instance, the order parame- 
ter critical exponent (3 takes a value 1/2) have not been 
convincingly verified numerically (a value [3 ~ 0.33 is 
systematically found in 2-dimensional Monte Carlo sim- 
ulations of the DLG and variations of it |12 |Jl3l), has 



triggered further studies. These new analyses ended up 
with the proposal of a new Langevin equation aimed at 
describing the critical properties of the DLG: 

5t./)(r,i) -TiVi0- Vi0+^Vi(/)3 + ri|V|0 + 7?. (2) 

This equation is a well known one: it coincides with 
the Langevin equation representing the random DLG 
(RDLG) (for which the driving field takes values 
00 and —00 in a random unbiased fashion, generating 
anisotropy but not an overall current). This equation 
has been extensively studied in flj]; its critical dimen- 
sion is dc = 3 (instead, dc = 5 for (|l|)) and its associated 
critical exponents and finite size scaling properties are 
now well known. At least two different recent numerical 
studies show that this equation reproduces rather nicely 
all DLG critical properties, and support the conclusion 
that it is the anisotropy and not the overall current the 
main ingredient characterizing the DLG phase transition. 

However, the situation is far from satisfactory. The 
central issue is that naive power counting analysis shows 
that the current term establishing the difference between 
the two abovementioned theoretical descriptions is a rele- 
vant perturbation at the Eq. (^ RG fixed point and there- 
fore, it is unclear why it should be absent in the proper 
Langevin description. It has been argued in ||l5|,|l2| that 
the coefficient of this term happens to vanish in the limit 
E ^ 00. This would imply that the DLG has a sort of 
multicritical point in the infinite fast driving limit. This 
scenario still needs to be confirmed numerically | |l^ , p^ . 

Another theoretical way out of this puzzling situa- 
tion is that the non-linear current term should be absent 
due to the fact that the microscopic theory is fermionic 
(i.e. occupation number restricted to be or 1), while 
the Langevin equation describes, in principle, a bosonic 
process: the current term is required in this bosonic 
formalism in order to have a vanishing current in per- 
fectly ordered aligned-with-the-field stripes (for which, 
the fermionic restriction precludes the presence of a non- 
vanishing current). We are presently working in the 
derivation of a field theoretical description that takes 
properly into account the fermionic nature of the micro- 
scopic model. 



III. SYSTEMS WITH MANY ABSORBING 
STATES 



Maybe the best well-known genuine non-equilibrium 
Langevin equation is the, so-called, Reggeon field theory 

dtpix, t) = Vp(x, t) + ap(x, t) - 6p^(x, t) + VP'?(x, t) 



(3) 



that captures the critical properties of phase transitions 
into a single absorbing state (with no extra symmetries 
nor conservation laws), usually referred to as the directed 
percolation (DP) universality class |^jr^;|l^. The key 
property of Eq. (||) is that all terms (including the noise) 
vanish in the absence of activity, i.e. for p(x) = 0). Even 
though convincing experimental realizations of this broad 
universality class are still missing, an overwhelming num- 
ber of models have been studied, all of them sharing their 
critical behavior with this minimal Langevin equation. 
The situation is less satisfactory for systems with many 
different absorbing states |l9), a prototype of which is 
the Pair Contact Process (PCP) [|0|j2l). In the PCP 
pairs of particles can generate new particles or get anni- 
hilated, but isolated particles do not have any dynamics; 
any configuration with just isolated particles is therefore 
absorbing. Models of this are relevant, for example, in 
catalysis The following Langevin equation for the 

PCP and related models was proposed some years back: 

5tp(x, t) = D\/^p{x, t) + ap{x, t) - bp^{x, t) 

+ ap{x,t)e~-^I>^-^^'^'' + ^prjix,t) (4) 

where the field p in Eq. ^ represents the density of pairs 
(activity), and the effect of the isolated particles (char- 
acterizing the different absorbing states |2^]) is captured 
in the non-Markovian exponential term. It has been ar- 
gued that the critical properties of this equation when 
approaching the critical point from the active phase are 
DP like | |2l| , p^ . Indeed, it is straightforward to see that 
in the presence of non-vanishing stationary activity the 
exponential term cancels out and we are left simply with 
DP. On the contrary, for spreading experiments for which 
critical propagation can occur inside the absorbing phase, 
the exponential term can be expanded in power series, 
and one ends up with |2^ 



dtpix, t) = D2Vlp{x, t) + ap{x, t) 
+ t) /q p(x, s)ds + ^77(x, t) 



(5) 



which is the well-known Langevin equation describing 
isotropic percolation dynamically i.e. dynamical percola- 
tion (DyP) These results are rather satisfactorily 
reproduced in numerical simulations pl|. Still there is a 
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point which remains obscure: If one works right at the 
critical point of the full theory, dynamical percolation 
terms are generated perturbatively, as first observed in 
pl| and recently stressed in |2j]. This new vertex be- 
ing more relevant than the dominant non-linearity in Eq. 
(^, leaded van Wijland to propose that the true asymp- 
totic critical behavior should be controlled by a DyP fixed 
point. In order to generate an active phase (missing in 
DyP) he proposes to treat the term proportional to in? as 
a dangerously irrelevant operator, and finds an analytical 
expression for j3. We believe that such a calculation can- 
not apply to the PGP since, even including the new term 
lacks of a well defined active phase. Being more precise, 
a term —hp cannot compensate the linear in time diver- 
gence of ap{x., t) /3(x, s)ds in the active phase. 

Another open problem in this context is understanding 
within a systematic RG calculation how the background 
field (describing the different absorbing configurations) 
emerges as a slave mode of the activity field, i.e. how it 
inherits the critical properties of the order parameter in 
the active phase p^ . A comprehensive understanding 
of this family of phase transitions putting together the 
active and inactive phases is still missing. 



IV. SELF-ORGANIZED CRITICALITY 

The observation that sandpiles, the archetype of self- 
organized systems |^ , fall into different absorbing states 
after every avalanche, right before new sand is added, 
opened the door to rationalize their critical properties us- 
ing Langevin equations similar to those described in the 
preceeding section. The first step in order to do so was 
to regularize the sandpiles, by introducing the, so-called, 
fixed energy sandpiles (FES) which eliminate dissipation 
and addition of energy (sand-grains) |p6|| . This converts 
the total energy into a control parameter: large total en- 
ergy generates stationary activity, while small amounts 
of energy lead the system with certainty to an absorbing 
configuration. The proposed set of equations for FES 
are: 

atp(x, t) = ap{x, t) - bp^ (x, t) + V2p(x, t) 

+ wE{x, t)p{x, t) + ^/pr^{:i^, t) 
dtE{K,t)^\\/^p{x,t) 



(6) 



where a, b, w, and A are constants, and r/ is a Gaussian 
white noise. In these equations the activity dynamics is 
controlled by the same type of terms appearing in Eq. (^ , 
plus an additional coupling between the activity field and 
a static conserved energy field. This extra term stems 
from the fact that creation of activity is locally fostered 
by the presence of a high background field density, and 
the energy is a conserved field, -^(x, t). The extra conser- 
vation law is therefore a new (relevant) ingredient with 
respect to RET. Some other terms, consistent with the 
symmetries and conservation laws, could have been in- 
cluded in Eq.(O) but they all turn out to be irrelevant 



from a power counting analysis |2^. This same set of 
Langevin equations (plus higher order noise terms) has 
been derived using Eock-space techniques for other dis- 
crete models with many absorbing states and a static 
local conservation law . The field theoretical analysis 
of this set of equations turns out to be a delicate issue 
(observe that analogous field theories for models with ab- 
sorbing states but where the conserved field is not a static 
one can be studied perturbatively without any problem, 
[p8|). As happens in the case of many absorbing states 
(without a conservation law) also here, at criticality DyP 
type of terms are generated. Here, even the physics com- 
ing from the active phase is not easy to work out. In 
this context, it has also been recently proposed Q that 
the critical properties should be described by the "reg- 
ularized" DyP fixed point, and again similar criticisms 
as those made before could apply here (although in this 
case, the problem is even more involved). 

A successful RG calculation of this theory would be ex- 
tremely valuable from a theoretical perspective; it would 
not only determine the critical exponents for a vast class 
of self- organized systems, but also clarify the issue of 
the proposed connection between self-organized critical- 
ity and the pinning of interfaces in disordered media [ p9[ . 

Another related problem is that of the study of the 
effect of quenched disorder in systems with absorbing 
states A field theoretical analysis by Janssen [^ , 

revealed the existence of running away RG trajectories, 
whose correspondence with the observed phenomenology 
in d = 1 and d — 2 |3l| | remains mysterious. 

Before finishing this section, we want to point out that 
a promising formulation of the same problem, namely de- 
riving an effective action for sandpiles has been recently 
addressed in . 



V. OTHER REACTION DIFFUSION SYSTEMS 

In this section we briefly enumerate some other open 
problems in field theoretical analyses of general reaction 
diffusion processes. 

A. Two symmetric absorbing states 

Eor some time it was believed that parity conservation 
(PG) was the main ingredient of a new, non DP, uni- 
versality class By now, it is well established that 
the presence of an exact Z2 symmetry between to equiv- 
alent absorbing states is its main distinctive trait 
Also, the introduction of parity conservation has been 
shown to play no relevant role in reaction-diffusion bi- 
nary spreading models ( 1 34 1 , see also ||3^] ) . A field theo- 
retical description of this universality class was proposed 
by Gardy and Tauber some years back. It starts with a 
Eock-space representation of the reaction diffusion lattice 
model in this class: A + A ^ 0, A ^ {m + 1)A (with 
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m an even constant) allowing to derive a field theoretical 
action. Even though the analysis of such a theory (that 
guarantees that the parity in the number of particles is 
conserved) is based on some uncontrolled expansion, it 
reproduces nicely many general features of this family, 
including the existence of a non-trivial critical point be- 
low two dimensions, the critical dimension dc — 2, and 
some other properties. The formalism is also applicable 
to odd values of m where it also generates sound results 

®- . . 

An interesting open problem in this perspective would 

be to construct a more adequate field theoretical repre- 
sentation that should include as the main ingredient the 
presence of two symmetric Z2 absorbing states, in the 
hope that in this more natural language a detailed stan- 
dard RG procedure would be applicable. 

Another related problem is that of writin g d own and 
analyzing a field theory for the Voter model and ex- 
tensions of it, for which a similar symmetry between dif- 
ferent absorbing states appears (also the non-equilibrium 
kinetic Ising model at zero temperature belongs two this 
family of models with Z2-symmetric absorbing states 

il). 



B. Pair Contact Process with diffusion (PCPD) 

In recent years, a new single-component absorbing- 
state universality class has been unveiled. It is the so 
called PCPD: if in the standard PCP we allow for dif- 
fusion of isolated particles we are led to this new class. 
Observe that switching-on diffusion represents a singular 
perturbation as, for instance, the many PCP absorbing 
configurations are reduced just to two, an empty one and 
one with a single wandering particle. It seems that the 
main ingredient in this class is the fact that reactions are 
binary (two particles are required for reactions to occur) 
and solitary particles travel performing random walks in 
between reaction zones. A field theory for this model 
was worked out by Howard and Tauber | |3^ some years 
back. Using a bosonic field theory (exact for a version 
of the model without a stationary active phase, usually 
called annihilation-fission process) they concluded that 
the critical dimension is dc = 2 and that the transition 
is not DP-likc. Unfortunately, the theory turned out to 
be non-renormalizable (i.e. an infinite hierarchy of rel- 
evant operators are generated perturbatively, making it 
un-tractable) . 

Different proposals have been made recently in order 
to rationalize the critical behavior of this class. These go 
from the existence of continuously variant exponents (as 
a function of the diffusion constant), to the existence of 
two universality classes (for small and large diffusion con- 
stants respectively), or just one well-defined set of critical 
exponents [ p9[ . 

One possible strategy to analyze this problem from 
a field theoretical point of view is to introduce dis- 



crete models in this class with two different species: one 
corresponding to diffusing "isolated" particles, and one 
"diffusing-reacting" type of particleplaying the role of 
the pairs in the original model |^,^. This leads to a 
set of Langcvin equations analogous to those proposed 
for the PCP, Eq.(|j) but including diffusion of the sec- 
ondary field. This changes the critical dimension from 
dc = 4 to dc = 2, but a systematic perturbative analysis 
allowing for a determination of the critical exponents has 
not been completed so far. 

Similar problems are observed upon studying ternary- 
reactions (as 3 A — > 0, combined with 3 A — > 3 + m) 
for which new critical behavior is expected For 
higher order nth-reactions the upper critical dimension 
for annihilation is below d = 1 therefore no anomalous 
phase transition is expected to occur [^J35| . 



VI. DISCUSSION 

We have briefly reviewed some non-equilibrium field 
theoretical open problems. They are related to non- 
equilibrium Ising-like models as well as systems with ab- 
sorbing states. Other families of problems not discussed 
here are interfacial growth, non-equilibrium wetting phe- 
nomena, transitions described by the multiplicative noise 
equation as for example those occurring in the synchro- 
nization of coupled-map-lattices, etc. The existence of 
the various problems reported here gives raise to the fol- 
lowing priorities for the developing of a systematic non- 
equilibrium field theoretical formalism: 

i) the necessity of developing new tools for deriving 
Langevin equations (or field theories) in a systematic, 
rigorous way, from discrete microscopic models. 

ii) Understanding the role of hard-core repulsion 
and/or implementing this constraint in a systematic way 
in field theoretical descriptions ji^ . 

iii) Developing new analytical schemes, specially in low 
dimensions, to deal with problems for which standard 
epsilon-expansion does not yield satisfactory results. 

It is my hope that this brief overview will stimulate 
further studies in this field. 
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